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Abstract. In this paper we consider the hyperbolic-elliptic system of two conservation laws 
that describes the dynamics of an elastic material governed by a non-monotone strain-stress 
function. Following Abeyaratne and Knowles, we propose a notion of admissible weak solution 
for this system in the class of functions with bounded variation. The formulation includes an 
entropy inequality, a kinetic relation (imposed along any subsonic phase boundary), and an 
nucleation criterion (for the appearance of new phase boundaries). We prove the ^-continuous 
dependence of the solution to the Riemann problem. Our main result yields the existence and 
the stability of propagating phase boundaries. The proof is based on Glimm's scheme and in 
particular on techniques going back to Glimm-Lax. In order to deal with the kinetic relation, 
we prove a result of pointwise convergence for the phase boundary. 

0. INTRODUCTION 

This paper deals with the following system of two conservation laws which describes the 
motion of an elastic material 

(0.1) d t w — d x v = 0, d t v — d x a(w) = 0. 

Here w > — 1 and v represent the deformation gradient and the velocity of the material, respec- 
tively. The stress function a :] — l,oo[— ► R is assumed to be monotonically increasing except 
in an interval ]wM,w m [. Such a form of the stress function is typical in the modeling of solid 
materials which admit different phases. A van der Waals gas also is described by a very similar 
system. System (0.1) is of mixed type, i.e. hyperbolic in the phase 1 region TLi = {w < wm} arid 
in the phase 3 region H.3 = {w > w m }, but elliptic in the intermediate region of phase 2 states. 
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The phase 2 states are known to be both mathematically and physically unstable (James [23]). 
We will consider here exclusively solutions which take their values in the phase 1 or phase 3 
regions only. 

Solutions to (0.1) in general are discontinuous and so must be understood in the sense of 
distributions; see Lax [25], [26] for background on weak solutions. Such discontinuous solutions 
are in general non-unique and those having a physical meaning must be selected through an 
admissibility (or entropy) criterion. We refer to Dafermos [7] for a review of entropy conditions 
in the setting of hyperbolic problems. As was pointed out by James [23], the mixed system 
(0.1) possesses a high degree of non-uniqueness, that a number of authors have attempted to 
resolve by means of suitable generalizations of entropy criteria from the theory of hyperbolic 
conservation laws. First of all, Shearer [38] considered the Lax entropy criterion [25], [26]. The 
viscosity and viscosity-capillarity approches have been analyzed by Slemrod [41], [42]; cf. also 
Hagan-Slemrod [18], Pego [35] and Shearer [39], [40]. Next, Hattori [19], [20] has investigated 
the application to (0.1) of the entropy rate admissibility criterion proposed by Dafermos [6]. 
Hsiao [22] has considered the Liu entropy criterion [32] which allows one to treat equations of 
state losing genuine nonlinearity in hyperbolic regions. Another approach to resolve the non- 
uniqueness can be found in a work by Keyfitz [24]. Additional material on system (0.1) is found 



All the above works consider the Riemann problem only, i.e. a Cauchy problem for (0.1) with 
initial condition which consists of two constant states. This problem can be solved explicitly 
(in a possibly non-unique way) by using simple waves (i.e. shock waves, rarefaction waves or 
contact discontinuities). Adding an "admissibility criterion" allows one to reduce the class of 
(admissible) solutions and in most situations to select a unique solution. However, it must be 
emphasized that the solution of the Riemann problem (when it is unique) depends on the chosen 
admissibility criterion. It turns out that there is no preferred criterion for the selection of the 
"physically meaningful" solutions of (0.1). 

A different approach was recently investigated by Abeyratne and Knowles in [2]. The main 
suggestion of these authors is that system (0.1) is not physically complete enough to describe 
the evolution of a phase boundary in an elastic material. It must be completed with a kinetic 
relation imposed along any subsonic phase boundary: this kinetic relation actually yields the 
rate of entropy dissipation across the phase discontinuity. Moreover, Abeyaratne and Knowles 
add an initiation criterion which controls the possible appearance of a new phase. We refer 
to [1] and the references therein for the motivation of introducing a kinetic relation and an 
initiation criterion which are actually classical in the context of quasi-static problems. Cf. also 
Gurtin [17] and Truskinovsky [44] for related ideas. 

Abeyaratne and Knowles proved in [2] that the Riemann problem for (0.1) always admits a 
unique admissible solution, i.e. a weak solution satisfying the kinetic relation and the initiation 
criterion, as well as the entropy inequality which reads 



in [12], [13] and [36]. 



(0.2) 



d t {W(w) + —) - d x (a(w)v) <0 



where W : ] — 1, oof— > R is the internal energy function defined by 




•w 



(0.3) 



a(y)dy for we] — 1, oof. 



Next they showed in [3] that the solution of the Riemann problem found by Slemrod through 
the viscosity-capillarity approximation corresponds to a special choice of kinetic relation in their 
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approach. It is not difficult to check also that the solution found by Shearer using Lax entropy 
inequalities coincides with the maximally dissipative kinetic relation investigated in [4]. (I thank 
Michael Shearer for pointing that out to me.) 

The present paper is devoted to continuing the analysis of system (0.1) through the approach 
of Abeyaratne and Knowles. As in [2], we will restrict ourselves to the case of a piecewise linear 
stress-function. This assumption simplifies the calculations but it is not a real restriction to the 
results of this paper. 

Our purpose is first (Sections 1 and 2) to give a slightly different presentation of the ideas of 
[2] , which as we think clarifies the concepts of kinetic relation and initiation criterion introduced 
by Abeyaratne and Knowles. Section 1 presents the mathematical formulation of a well-posed 
(at least for Riemann data) problem associated with system (0.1). As is usual for hyperbolic 
problems, we consider bounded solutions of bounded variation (BV). Our formulation follows 
[2] with however two main modifications. The kinetic relation is introduced from a completely 
dynamical point of view and not as a generalization of the quasi-static point of view as was 
done in [2]. That leads us to a larger range of admissible values for the — as we call it below 
- entropy dissipation function in the kinetic relation. Furthermore, the initiation criterion at 
some point x is formulated here in two different ways, depending on whether x is in an interior 
point of the space interval [a, b] where we set the problem, or x is a point of its boundary. For 
definiteness, we allow spontaneous initiation of a new phase only at the extremities of the bar 
[a,b], which is consistent with the classical static theory. 

Then Section 2 describes briefly the solution of the Riemann problem. We explain how to 
take into account the two changes above in the construction of [2]. The main result of this 
section establishes the L 1 continuous dependence of the Riemann solution with respect to its 
initial states. It must be emphasized that the two observations above are essential for the 
continuous dependence property to hold, especially our condition that a new phase may occur 
spontaneously only at the extremities of the bar. The same results are also obtained for the 
Riemann problem in a half-space. Note that, although uniqueness of the admissible solution 
holds for the Riemann problem, nothing is known for the general Cauchy problem. As a matter 
of fact, the issue of uniqueness for conservation laws is understood in a few number of situations 
only. (See, for hyperbolic problems, LeFloch-Xin [31] and the references therein.) 

The second part of the paper (Sections 3 and 4) focuses on the solutions of the Cauchy 
problem for system (0.1), which are BV perturbations of a single propagating phase boundary 
separating a phase 1 state and a phase 3 state. We prove the existence of admissible weak 
solutions of this form, when the initial data on both sides of the phase discontinuity has small 
total variation. We treat the case of any non-characteristic phase boundary as well as the case 
of a characteristic phase boundary provided that no strong wave arises from perturbating the 
states on both sides of the phase boundary. The random-choice scheme due to Glimm [15] is 
used to construct approximate solutions to the problem. Its stability in the BV norm is proved 
from an essentially linear estimate of wave interactions between two Riemann solutions. Such 
linear interaction terms were used in a different situation by Chern [5] and Schochet [37]. Note 
that the strength of the phase discontinuity is not (and can not be) assumed to be small in any 
sense. 

The stability of the scheme in the total variation norm is sufficient to extract a subsequence 
converging to a weak solution of the problem. This convergence result holds almost everywhere 
with respect to the Lebesgue measure. This is sufficient to show that the scheme converges to a 
weak solution of the problem. But, proving that this solution is admissible requires a result of 
pointwise convergence of the phase boundary. In Section 4, we establish this property by using 
the technique of analysis due to Glimm-Lax [16]. We next prove that it is sufficient, at least for 
non-stationary phase boundaries, for the passage to the limit in the kinetic relation. 
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An extension of the results in this paper to arbitrary large initial data would require a better 
understanding of the phenomena of initiation of new phases. 

Many ideas in this paper are related to those in the developing theory of nonlinear hyperbolic 
systems in non-conservative form for which we refer the reader to Dal Maso-LeFloch-Murat [8] 
and LeFloch-Liu [30]; see also [27] to [29]. 
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1. MATHEMATICAL FORMULATION OF THE PROBLEM 

This section describes the formulation of the Cauchy problem associated with the mixed 
system (0.1). The formulation includes the system of conservation laws (mass, momentum) (0.1) 

together with the (Clausius-Duhem) entropy inequality associated with the entropy W(w) + \- 
It is made complete by adding to these both a kinetic relation along any subsonic phase boundary 
and an initiation criterion for the occurrence of possible new phase boundaries in the solution. 
We specify below the assumptions on the kinetic relation and the initiation criterion which will 
be essential to the results of Section 2. This section also introduces notation which will be of 
constant use throughout this paper. 
We write system (0.1) in the form 

(1.1) d t u + d x f(u) = 0, u=( V ), f(u)- ( ~ n[ ' n 



w J \ —v 

For simplicity, we shall assume that the stress- function a:] — l,oo[— > R is a piecewise linear 
function of the following form 

k\w for — 1 < w < wm, 

(1.2) <t{w) = { k 3 w m + (hwM ~ hw m )(w - w m )/(w M ~ w m ) for w M <w <w m , 
k 3 w for w m < w. 

The constants fei, ks,w m and wm in (1-2) are assumed to satisfy the properties 

(1.3) < < ki and < wm < w m . 
We shall use the notation 

a M = k x w M and a m = k 3 w m . 

The phase 1 region 7i\ = { — 1 < w < u>m} and the phase 3 region TC 3 = {w > w m } correspond 
to observable and stable states. In our formulation below, the solution cannot enter the unstable 
phase 2 region {u>m < w < w m } and so must jump from Hi toHs or conversely. A discontinuity 
between two states in different phases is called a phase boundary. 
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System (1.1) is linear hyperbolic in Hi and H3, and the corresponding characteristic speeds 
are 

±Ci = ±y/ki in Hi and ±03 = ±\/fc3 in H3. 

In view of (1.3), the waves in phase 1 travel faster than those in phase 3, i.e. c\ > C3. We may 
also use the notation 

(1.4) c(w) = C\ if w < wm, C3 if w > w m . 

Note that c(w) is not defined if w belongs to ]wm ,w m [. With some abuse of notation, Hi and 
H3 will sometimes also denote {(«, w)/ — 1 < w < Wm, v G R} and {(v,w)/w m < w, v G R}, 
respectively. We also set H = Hi U H3. 

Since (1.1) is linear hyperbolic in H\ and H3, possible discontinuities in the initial data for 
(1.1) are simply advected along the characteristic lines of slopes either ±ci or ±03. (This is 
true at least up to the time of appearance of a new phase.) The elementary waves in each 
of the regions Hi and H3 are contact discontinuities. Hence, the special choice (1.2) for the 
constitutive law is very convenient. It makes quite simple the analysis in the hyperbolic regions 
and allows us to focus on the phase boundaries between Hi and H3. We will see that the 
description of the appearance and the evolution of the phase boundaries is far from trivial. 

In the theory of hyperbolic conservation laws, it is standard to consider solutions u = (v, w) 
to (1.1) in the functional space L^ c (R + x R, H) (recall that H = Hi U H3) which satisfy 

(1.5a) system (1.1) in the sense of distributions, 



(1.56) the entropy inequality (0.2), (0.3) in the sense of distributions 

and 

(1.5c) an initial condition -u at t = in the L\ oc sense. 

Here no is a given function in Lf£ c (R,H) and, for future reference, we rewrite the entropy 
inequality in the form 

(1.6a) d t U{u) + d x F(u) < 

with 

v 2 r 

(1.66) U(u) = W(w) + — , F(u) = -a(w)v and W(w) = <r(y)dy. 

2 Jo 

We recall that solutions in the sense (1.5) are unique (at least for Riemann data) in the standard 
situation of a (genuinely nonlinear or linearly degenerate) increasing strain-stress function a. 
This is no longer true in the case of the mixed system under consideration here: see for instance 
James [23]. We also point out that the entropy function U is not a convex function. 

To complete the formulation (1.5), we follow Abeyaratne and Knowles in [2]. Let us give 
first some motivation for their suggestion. Suppose u e = (v e ,w e ) is the solution of a regularized 
version of system (1.1) obtained by adding high-order terms, depending on a (small) parameter 
e, in the right-hand side of the equations (e.g. use the viscosity-capillarity terms as was done 
by Slemrod [41]). As was pointed out by Lax for general systems of conservation laws, the limit 
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u = lim u e - if it exists (and if the convergence holds in a suitable topology) - must be a solution 
to (1.1) in the sense (1.5); in particular the entropy inequality (1.6) must hold. Since (1.5) is 
incomplete, it seems natural to "keep more information" about the limiting function u from its 
regularization u e . Specifically Abeyaratne and Knowles' suggestion is equivalent to replacing 

(1.6) with the stronger requirement that 

(1.7) dtU(u) + d x F(u) = ti, 

where \i is a given non-positive measure that clearly must satisfy certain restrictions. Note that 
in principle jj, could be determined by the formula 

H = weak-star lim(d t U(u e ) + d x F(u e ) ) 

(at least when u e has uniformly bounded total variation in (t, x)). This formula may not give a 
very explicit expression for fi. Hopefully, it turns out that (1.7) is needed (to achieve uniqueness) 
only for one kind of discontinuity: the subsonic phase boundaries. Moreover, in that case, we 
can allow a large range of measures \i. Here, we call subsonic (respectively supersonic) those 
phase boundaries that travel with speed less (resp. greater) than the contact discontinuities in 
phase H3. 

The precise formulation of condition (1.7) given below requires that u is a bounded function 
of bounded variation. When u has bounded variation, we call entropy dissipation the value of the 
measure d t U(u) + d x F(u) along a curve of (contact or phase) discontinuity of u. According to 
[2], the kinetic relation yields this entropy dissipation along any subsonic phase boundary, as an 
explicit function, say (/>(V), of the speed V of propagation of this discontinuity. In applications, 
the actual kinetic relation, that is the function cf), must be determined from the properties of the 
specific material under consideration. This kind of constitutive model is already in extensive 
use in the quasi-static setting for problems of phase transition in solids. We refer the reader 
to [1] as well as Truskinovsky [44] and the references cited there. The speed V can also be 
interpreted as an internal variable and the kinetic relation indeed determines the evolution of 
this internal parameter. 



Remark 1.1. 1) That subsonic and supersonic phase boundaries must be treated in a different 
way is clear, for instance when solving Riemann problems. A wave structure with a supersonic 
phase boundary contains two waves, while one with a subsonic boundary is composed of three 
waves. This latter case suffers, without a kinetic relation, from a strong lack of uniqueness. 
Cf. James [23] and Section 2. 

2) The approach considered here has some similarity to the theory of nonlinear hyperbolic sys- 
tems in non-conservative form; cf. Dal Maso-LeFloch-Murat [8] and LeFloch-Liu [30]. Namely, 
as is the case for systems (1.1), the weak solutions to these systems are not uniquely determined 
by the partial differential equations and an entropy inequality, but an additional constitutive 
relation must be added to ensure uniqueness. This fact was first pointed out by LeFloch; cf. [27] 
to [29]. 

3) Conservation laws with measure source-term like (1.7) have been useful in various contexts, 
cf. Di Perna [9], Di Perna-Majda [11], Hou-LeFloch [21]. 

Let us introduce some notations and recall some facts about functions of bounded variation, 
that can be found in Volpert [45] and Federer [14]. Let Q be an open subset of R m . A func- 
tion u-.n->R p belongs to the space BV(Sl,R p ) (respectively BVi oc (n,R p )) Hue L 1 (Q,R P ) 

du 

(resp. L} oc (tt,R p )) and the distributional derivatives — — for 1 < j < m are bounded (resp. 
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locally bounded) Borel measures on Q. In what follows, we will always consider functions in 
L°°(n,R p ) n BV(n,R p ) or L% C {Q,R P ) n BVi oc {Q,R p ), often called for short BV functions or 
BVioc functions. For each BVi oc function u, we have the following decomposition 

n = c{u)us{u)uE(u), 

where 

C(u) is the set of all points of approximate continuity for u, 
S(u) is the set of all points of approximate jump for u 
and 

E{u) is the set of exceptional points with the property H m _i(E(u) ) = 0. 
Here i? m _i is the (m — 1) -dimensional Hausdorff measure on R m . For each point y in S(u), 
there exists a unit normal v 6 R m and approximate left and right limits for u that we denote 
by u±(y). The set S(u) consists of the union of a countable number of rectifiable curves. 

We denote the norm of u by |M|w(n,,R p ) = IMlL 1 (f2,,R p ) + \Du\({l), where Du is the measure 
du du du 

(— , — , - — ). When u = u(t,x) G L^(R + x R,H) D BV loc (R + x R,H), we use the 
dyi dy 2 dy m 
notation: 

v(t,x) = (y t (t,x), v x (t,x)) and V(t,x) = - Ut ^ x ^ 

u x (t,x) 

valid for all (t, x) £ S(u). The ratio V(t,x) represents the speed of propagation of the disconti- 
nuity in u at the point (t, x). Note that system (1.1) has the property of propagation with finite 
velocity (in regions TL\ and H3). So v x (t, x) will never vanish, and for definiteness we always 
choose u x (t,x) > 0. In the following, we shall always have: u(t) £ BV for all times t. 

Let (p: } — C3, cs[~ > R be a function, called below entropy dissipation function, satisfying the 
following properties: 

(1.8a) (p belongs to C 2 (] - c 3 ,0[U]0, c 3 ]) and 0(O±) and 0'(O±) exist, 

(1.86) lim = ip(cs) and <j>'"{cz—) exists, 

v^c- 

(1.8c) lim = —00, 

(1.8d) is increasing on ] — C3, C3] 

and 

'^(y)<0(y)<o for fg]-c 3 ,o], 



(1.8e) 



< <f>{V) < ip(V) for V € [0, c 3 ]. 



In (1.8b) and (1.8e), the minimal and maximal entropy dissipation functions ip: ] — c 3 , 0] — > i?_ 
and ip: [0, c 3 ] — > R + are defined by 

(1.9a) ±{V) = {kl ~ fc3 W (^ ~ ^ I for ^e]-c 3 ,0] 
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and 
(1.96) 



V{V) = w m (w M 



h-V 2 

h-v 2 



w m ) for V£[0,c 3 }. 



Remark 1.2. 1) Inequalities (1.8e) give the range of values taken by the entropy dissipation 
rate £{u) (see below) when varying the left and right values at a discontinuity satisfying the 
Rankine-Hugoniot relations and the entropy condition. 

2) In [2], instead of (1.8e), Abeyaratne and Knowles assume the (more restrictive) condition: 



(1.8e)' 



' ^{0)<(f)(V)<0 for Ve]-c 3 ,0], 



< (f>(V) < -0(0) for Ve[0,c 3 ]. 



3) Assumptions (1.8) made in this paper are indeed satisfied in the examples considered by [3] 
and [4]. For instance, they are fulfilled by the maximally dissipative function max defined by: 



.(V)=i>(V) for V€]-c 3 ,0], i/>(V) for Vg[0,c 3 ]. 



We next define the entropy dissipation rate £(u) associated with any function 

u € L% C (R+ xR,H)n BVi oc {R + x R,H) 

by the following formula 

(1.10) £(u) = -(U(u + ) - U(u.) ) - ^(F(u + ) - ) 

which defines £(u)(t,x) at i^i-almost every point (t,x), where v t {t,x) ^ (i.e. V(t,x) ^ 0). 
£(u) is the product of — ^- by the jump of the measure d t U(u) + d x F(u) along the curve of 
approximate jump of u. Formula (1-10) makes sense only if v t (t, x) 7^ 0. However, it is a simple 
observation that if u is assumed to be a weak solution to system (1.1), then the above jump (i.e. 
the entropy dissipation) vanishes at the points where v t vanishes. This fact allows us to define 
£(u)(t,x) Hi -almost everywhere, as shown by the following lemma. 

Lemma 1.1. Ifu G Lf^ c (R + x R,H) n BVi oc (R + x R,H) is a weak solution to (1.1), then one 
has 

(1.10)' £{u) = - f W+ |a(y) - \{a(w + ) + *(«;_))} dy, 

at Hi-almost every (t,x) such that f t (t, x) 7^ 0. 
From now on, we use (1.10)' to define £(u)(t, x). 
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Proof. At a point of approximate discontinuity (t, x) of the solution u, the following Rankine- 
Hugoniot relations hold: 

v t (w + - w-) - u x (v + - v-) = 0, 

v t (v + - v-) - u x (a(w + ) - a{w-) ) = 0. 

These relations used in (1.10) yield: 

r w + i Vx 

-£(u) = / a(y) dy + -(v + - v_) - —(a{w + )v + - a(w-)v-) 

r w + i Vx Vx 

= I v(y)dy + -(v + + v_) — (a(w + ) - a(w-)) - —(a(w + )v + - a(w-)v-). 

We thus get 

r w + Vx 

—£(u)= / a(y) dy + — —{v + a(w + ) + V-a(w + ) - v+a{w-) - V-cr(w-) 
— 2a(w + )v + + 2a(w-)v_}, 

so that 

r w + i Vx 

-£(u) = / a(y) dy - - — (v+ - v-)(a(w+) + a(w-) ), 

which, in view of the Rankine-Hugoniot relations above, gives the desired result (1.10)'. □ 

Let us denote by B su b(u) the set of all points of approximate discontinuity in a weak solution 
u that correspond to a subsonic phase boundary. This means: 



B su b(u) = <j (t,x) G S(u)\ either : u-(t,x) G Hi, u+(t,x) G H3 and \V\ < C3, 
or : U-(t,x) G Hz, u + (t,x) G Hi and \V\ < C3 >. 



In view of (1.6), the Borel measure d t U(u)+d x F(u) is globally non-positive. The kinetic relation 
now specifies the value itself (and not only the sign) of this measure along any subsonic phase 
boundary. In other words, for iii-almost all (t, x) G B su b(u), one must have 



(1.11) £(u)(t,x) 



-<j)(V(t,x)) if U-(t,x)€Hi, 
4>(-V(t,x)) if U -(t,x)eH 3 . 



Remark 1.3. As a matter of fact, the traveling waves obtained through the viscosity-capillarity 
regularization to system (1.1) converge to weak solutions of (1.1) that satisfy the kinetic relation 
(1.1) with a specific choice of function <j). This function can be determined explicitly and depends 
only on the viscosity and capillarity coefficients introduced in the regularization (cf. [3]). 
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Finally, we have to formulate the initiation criterion, which together with the above kinetic 
relation will allow us to rule out all non-physical solutions to our problem. Let ]a, b[ be a space 
interval in which we are going to set the problem, with a < b and possibly a = — oo and/or 
b = +00. The initiation criterion will reflect the following facts: 



From the mathematical point of view, condition (1.12) is essential: it ensures that sponta- 
neous initiation of a new phase inside ]a,b[ cannot occur from two nearby initial states in the 
same phase (cf. Section 2). This does not exclude the possibility (and it really happens) that 
an initial discontinuity with large jump gives rise to, for instance, a phase 1 state although the 
states on both sides of the initial discontinuity are in phase 3. However, by condition (1.12), a 
single constant state is always a (trivial) admissible solution. (This property was not satisfied 
in the construction of [2].) This is also essential to get the L 1 continuous dependence property 
for Riemann solutions, proved below in Section 2. 

Conditions (1.13) and (1.14) follow the quasi-static theory [1]. They allow "spontaneous 
nucleation" of a new phase only at the end points of [a, b\. Note that, more generally, we could 
as well allow nucleation at some arbitrary given points of [a, b}. Our actual restriction is that 
the points of spontaneous nucleation are known a priori and follow a selection criterion of the 
form specified below. However, while this formulation is fully satisfactory from the mathematical 
point of view, it does not reproduce what is really observed in practical experiments with elastic 
bars. Namely, in experiments, when pulling out an elastic bar uniformly in phase 1, initiation 
of phase 3 regions in the bar occurs successively and (apparently) randomly at various places in 
the bar. Physicists assert that initiation occurs at microscopic inhomogeneities of the material. 
A complete treatement of the initiation mecanism is beyond the scope of this paper and would 
probably require a statistical description. (As a matter of fact, this might be included in the 
random choice scheme, studied below, quite easily.) 

It remains to provide an analytic version of the conditions (1.12) to (1.14). For convenience, 
we use here an averaged strain in our formulation. (In [1] and [2], the stress and the entropy 
dissipation rate, respectively, are used instead.) Given any function u = (v,w) in L^ c (R + x 
R, H) n BV loc (R+ x R,H), we set 



(1.12) 



no new phase occurs from any point x in ]a, b[ 

except if no solution exists without creation of a new phase, 



(1.13) 




and 



(1.14) 




(1.15) 



K = 



c{w-)w- + c{w + )w + + v + — V 



c(w-) + c(w + ) 
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which defines h u (t,x) for Hi - almost every (t, x) in R + x R. We note that h u (t,x) = w(t,x) 
when (t,x) is a point of approximate continuity of u. So h u (t,x) represents an averaged strain 
at the point (t, x) and determine the dynamics at this point. For instance, if u_ and u + are in 
the same phase, h u is the intermediate value between the 1-wave and the 2-wave in the solution 
of the Riemann problem with initial data u_ and u_|_ (cf. Section 2). 

At any interior point x G]a, b[ and for each time t > 0, the initiation criterion, by definition, 

is 



(1.16) 



' if U-(t,x) and «+(i, x) belong to TCi (respectively H 3 ), then: 

h u (t,x) G Hi (resp. H3) if and only if there exists e > such that 
k u(s, y) G Hi (resp. H3) for (s, y) e[t,t + e[x]x - e, x + e[. 



According to (1.12), condition (1.16) ensures that, locally in time, the solution remains in the 
same phase whenever this is possible. Cf. Section 2. 

We are now concerned with the boundary points x = a and x = b. We assume that u(t) is 
defined for all times and has bounded variation in x. (This will be the regularity of the solutions 
found in Section 4.) The material is assumed to be fixed at the end points, i.e. when a ^ —00 
and/or b 7^ +00, we have 

(1.17a) v+(t, a) = for L 1 - almost every t > 0, 

and 

(1.176) v-(t, b) = for L 1 - almost every t > 0, 

where L 1 denotes the one-dimensional Lebesgue measure. Since v has bounded variation, it 
admits a L 1 trace at x = a and x = b. Let w c ^ and be two constants, called critical values 
for the initiation, that must satisfy the inequalities 

(1.18) ^<wZ<w M and w m < wZ < ^, 

Kl K3 

where gq is the so-called Maxwell stress given by 

(1.19) cr = ^a m a M = cic 3 ^w m w M - 

Note that, as pointed out to us by Abeyaratne, the critical values for intitiation should in 
principle depend on the speed of propagation of the phase discontinuity. At the point x = a 
(when a 7^ —00), we impose for all times t > the following two conditions: 

{if u + (t, a) belongs to H3, then: 
h u (t, a) > wZ if and only if there exists e > such that 
u+ (s, a) G H3 for s G [t, t + e[, 

and 

if u + (t, a) belongs to Hi then: 

(1.20) jj ^ h u (t, a) G Hi if and only if there exists e > such that 
u + (s, t) G Hi for all s G [t,t + e[. 
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In order to satisfy the boundary condition (1.17a), the term h u (a,t) in (1.20)j is defined by 
formula (1.15) with here 

(1.21) V-(t, a) = — v + (t, a) and w_(t, a) = w + (t, a). 

Similarly, at the point x = b (when b ^ +oo), we impose for all times t > that 

{if U-(t, b) belongs to Hi, then: 
h u (t, b) < w°m if and only if there exists e > such that 
u- (s, b) G Hi for all s G [t, t + e[, 

and 

{if u-(t, b) belongs to Hz, then: 
h u (t, b) G H3 if and only if there exists e > such that 
u- (s, b) G H 3 for all s G [M + e[. 

As previously, we set here: 

(1.23) v+(t,b) = -v-(t,b) and iu + (t, b) = tu_(t, 6). 

We call an admissible weak solution to system (1.1) a function u = (v,w) which satisfies the 
conservation laws (1.1), the entropy inequality (1.6), the kinetic relation (1.11), the boundary 
condition (1.17) (if instead of R an interval ]a,b[ is considered) and the initiation criterion 
(1.16), (1.20) and (1.22). 

In this paper, we will prove the existence of such an admissible weak solution for two kinds of 
Cauchy data: the Riemann problem (in the whole space and in a half space) and a perturbation 
of a single propagating phase boundary. These results provide a strong justification for our 
formulation here. It would be interesting to address the general question of existence and 
uniqueness for system (1.1) in the setting introduced in this section. 

Remark 1.4. 1) A phase boundary necessarily is a wave with a large strength (at least 
\wm — w m \). This implies that, in BV solutions, phase boundaries cannot accumulate in a 
bounded region of the (t, x)-plane. Phase boundaries are thus isolated, and this justifies the 
formulations (1.16), (1.20) and (1.22). 

2) Sections 3 and 4 will provide an existence result for small BV perturbations of phase 
boundaries. We believe our result to be true for any finite number of phase boundaries. 
However, for arbitrary large data, the appearance of an infinite number of phase boundaries is 
not excluded a priori. In such a case, the solution would not have bounded total variation. A 
challenging issue is to extend the present formulation to the framework of solutions. 

3) The formulation of this section can be extended to the case that the stress-function is not a 
piecewise affine function but an arbitrary piecewise monotone function. 

4) Shearer's solution [38] corresponds to the choice w 1 ^ = Wm, = w m and 4> = </>max (see 
Remark 1.2 for the definition of ^ max )- 
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2. RIEMANN PROBLEM AND CONTINUOUS DEPENDENCE 

This section gives an explicit description of the admissible weak solution of the problem 
formulated in Section 1, in two cases: the Riemann problem in the whole space and the Riemann 
problem in a half space. Our main result in this section is the L 1 continuous dependence property 
of the solution of these problems. Note that the formulation of Section 1 and the assumptions 
made there are essential for this property to hold. 

We consider the following two problems: 

1) the Riemann problem in the whole space ]a, b[=] — oo, +oo[ which corresponds to an initial 
data of the form 

(2.1) uq(x) = ul for x < and ur for x > 0. 

Here ul = (vl,ujl) £ H- and ur = (vr,wr) £ TL are two given constant states. 

2) the Riemann problem in the half space ]a, b[=]0, +oo[ which corresponds to the initial data 

(2.2) Uq(x) = u for all x > 
where no = (vo,Wo) € TC is a constant state. 

We shall describe successively the admissible solutions to problems 1) and 2), by following 
closely the work by Abeyaratne-Knowles. However our construction is slightly different from 
the one in [2], due to our formulation. We will not address the question of uniqueness of the 
solution here, since it is an easy matter from the results in [2] (which yield for their construction 
uniqueness in the class of solutions composed of simple waves) . 

To begin with, we deal with problem 1) and distinguish between several cases: 

case l)-a: ul £ TCi and ur £ H3, 

case l)-b: ul € H-i and ur 6 Hi, 

case l)-c: ul H3 and ur 6 H3, 

case l)-d: ul G H3 and ur 6 H\. 
Cases l)-c and l)-d are very similar to cases l)-b and l)-a, respectively. (Use the transformation 
x — ► — x and the fact that the equations (1.1) and more generally all the requirements in the 
formulation of Section 1 are invariant under this transformation.) So we omit them and focus 
on the two first cases. 

Case l)-a: Suppose that ul 6 7i\ and ur 6 H3. 

We must construct a solution to (1.1), (2.1) which is admissible in the sense of Section 1. 
The solution necessarily contains a phase boundary (and only one as was checked in [2]) with 
phase 1 states at the left and phase 3 states at the right. Two different wave structures are 
possible, depending on whether the phase boundary is subsonic or supersonic. Let V be the 
speed of the phase boundary and set 

(2.3) h LR = h u (0,x) = ; {C!W L + c 3 wr +v R - v l ). 

C\ + c 3 

We distinguish between two cases depending on the sign of h^R. 

Case l)-al: Suppose moreover that h^R > 0. 

In that case, we seek the solution u in the following form 



(2.4) u(t,x) 





for 


X < —C\t, 


U- 


for 


- cit < x < Vt, 


< 

u + 


for 


Vt < X < Cj,t 


. ur 


for 


X > c^t, 
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where the constants U- = (v-,W-) and u + = (v + ,w + ) belong to Hi and H 3 respectively. The 
solution contains a contact discontinuity of speed — Ci, the phase boundary with subsonic speed 
\V\ < C3 and a contact discontinuity with speed c 3 . We are going to prove that indeed such a 
solution exists by determining explicitly the values of the constants U-, u + and V. 

First of all, u given by (2.4) must be a weak solution to (1.1), so must satisfy the following 
four Rankine-Hugoniot relations 



(2.5) 



{C\{w- — wl) — (v- — Vl) = 0, V(w + — W-) + v _|_ — V- = 0, 
V(v + - v-) + cjw + - c\w- = 0, c 3 (w + - w R ) + v + - V R = 0. 



If V is chosen as a parameter, (2.5) yields explicit expressions for V-, v + , W-, w + as functions 
ofV: 



(2.6) 



, C 3 + V h C 3 + V U 

v- = v L - ciw L H —c x h L R, w- = —h LR , 

ci + V ci + V 

v+ = v R + c 3 wr — c 3 h LR , w + = — h LR , 

c 3 -V c 3 -V 



Formulas (2.6) define a one-parameter family of solutions to problem (1.1), (2.1). Note that W- 
and w + are always non-negative. 

Next we take into account the kinetic relation that states (cf. (1.10) and (1.11)): 

U(u + ) - U(u_) + ^{F(u + ) - F(u.) ) = 0(F), 

or using the more general form (1.10)': 

" Ml/) - \H w +) + )}dy = 4>{V). 
Using the expression (1.2) for the function <r, this becomes 

(2.7) ^(fci - k z ) (w M w m - w + w-) = 4>{V). 

If we use in (2.7) the expressions for W+ and W- given by (2.6), it follows that 

(2.8) W = 0(n where 6(V) = ^^[w M w m - |g ± ^ = ^ 

Note that the function depends only on the averaged strain Hlr- In view of our set of assump- 
tions (1.8) and 

9' < 0, 0(-c 3 ) = ?(c 3 ) and lim 0(F) = -00, 

one easily checks that equation (2.8) admits a unique root V. Moreover, if this specific value of 
V is used in (2.6) to get u>_, w + , V- and v + , then our construction is consistent in the sense 
that one has 

(2.9) w- G Hi and w+ € Hz- 
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We now prove (2.9) by using the assumption (1.8e) on the function 4> (a stronger assumption 
was made in [2] to derive (2.9)). 

Let us for instance check that u>_ 6 Hi, in other words: W- < wm- In view of (1.8e) and 
(2.8), one has ip(V) < 6(V), so that by (1.9a) and (2.8) again 



ki-V 2 2 ^ (c 3 + V)( Cl -V) 

" " Wm < WMWm - (c 1 + V)(c 3 -V) " ! 



w M w m - ^ — w M < w M w m - )~° ' l^)" 1 h 2 

Since |V| < C3, we obtain 



w 2 >h 2 {Cl + V? 
But wm > and h LR > by assumption, so 

w M > h LR — = w- 

C3 + V 

in view of the expressions (2.6). The proof of (2.9) is complete. 

Finally, we note that the entropy inequality (1.6) is trivially satisfied along the contact waves, 
while it is a consequence of the kinetic relation (1.10), (1.11) along the phase boundary. Thus, 
in the present case, (1.6) yields no additional constraint. 

The above construction yields the admissible weak solution of the problem. Based on the 
explicit expressions (2.6) and the implicit equation (2.8), it is elementary to prove the following 
regularity result of the Riemann solution. 

Lemma 2.1. Consider the Riemann problem (1.1), (2.1) in case l)-al, i.e. with ul € Hi, 
u R 6 H3 and h^ R > 0. Then the admissible weak solution to this problem is given by formulas 
(2.4), (2.6) and (2.8). One can consider the states U- and u + and the speed V in (2.4) as 
functions of the initial states ul and u R , or more precisely 

v- = v- (u L , h LR ), w- = w- (h LR ) 

and 

v + = v + (u R ,h LR ), w + =w + {h LR ) and V = V(h LR ). 

The functions u_, u + and V are Lipschitz continuous in the range of values {ul 6 Hi, u R 6 
I Iilr > 0}. They are of class C 1 (with Lipschitz continuous derivatives) away from V = 0. 
The behavior of U-, u + and V when /i_l_r — > + is given as follows: 

(2.10a) lim v-=vl — c\Wl, lim u>_ = 0, 



(2.106) lim v+ = VR+Ca i WR -j£M\ lim w + - ■ ^ 



h L R^o+ 1 y c 3 j h LR ^o+ y C3 

(2.10c) lim V = c 3 , lim = (c 3 - a) 



h L R^0+ " h LR ^o+ dh LR V <j)'(c 3 

and 

(2.1M) lim #^ = ^i^, lim dw ~ - 



h LR ^o+ dh LR C1 + C3 h LR ^o+ dh LR C1 + C3' 
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Remark 2.1. 1) Assumptions (1.8a) and (1.8d) imply that 4> 1 exists and is a Lipschitz 
continuous function. Away from V = 0, _1 is of class C 2 , and so is the function V(1ilr) in 
view of (2.8). 

2) If the function (f> 6 C 2 (] — 03,03]), then all the functions in Lemma 2.1 are globally of class 
C 2 . (This is not the case of the maximally dissipative function quoted in Remark 1.2.) 

3) In the special case that 4> = <^max, we find: J - ^ = w m . 



UL 


for 


X < —Clt, 


U- 


for 


- cit < x <Vt 


ur 


for 


x > Vt. 



Case l)-a2: Suppose now that Hlr < 0. 

In that case, the solution is composed of a contact discontinuity with speed —c\ and a phase 
boundary with supersonic speed V > c 3 . There is no C3 contact wave. We use the notation 



(2.11) u(t,x) 



The state u_ = (v-,W-) € Hi and the speed V must satisfy the following jump relations: 
ci(w- - w L ) - (v- - v L ) = 0, 

V(w R - wJ) + v R - v- = 0, V(v R - v_) + cj-WR - c\w- = 
We thus get V- and W- explicitly as functions of V 

(2.12a) V- =Vl — C\Wl H — - cih LR , w- = — h LR , 

ci + V ci + V 

the speed V being given by an implicit algebraic equation 

(2.126) V 2 { - c 3 w R + (ci + c 3 )/i L i?} + (c| - c?)iWflV + (ci + 03)^(03^ - ci/i Lfl ) = 0. 

Note that w_ given by (2.12a) is always non-positive. 

One can check [2] that (2.12b) has a unique solution which belongs indeed to the (physically 
interesting) interval [03, Ci[ if and only if Hlr satisfies the restriction 

(2.13) h^KhLRKO with h 00 = — - — {c 3 w R -Ci-(clw 2 R + (ci + c 2 3 )wR + ci) 1/2 }. 

Cl + c 3 

In other words, the Riemann problem can be solved when Hlr < if and only if Hlr > h^. We 
emphasize that the kinetic relation was not used here and V is found to be supersonic; this is 
in complete agreement with the fact that (1.11) is imposed only for subsonic phase boundaries. 

Lemma 2.2. Consider the Riemann problem (1.1), (2.1) in the case l)-a2 that ul € TCi, 
ur € 7i 3 and h^R < (with the restriction (2.13)). Then the admissible weak solution to this 
problem is given by (2.11), (2.12). One can consider the state U- and the speed V as functions 
of the initial states ul and Ur, or more precisely: 

v- = V-(u L ,h LR ,w R ), w- = W-(h LR ,w R ) and V = V(h LR ,w R ). 
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Then the functions u_ and V are C°° functions of their arguments and when h^ R — > 0—, they 
satisfy: 

(2.14a) lim v-=vl — c\Wl, lim w- = 0, 

flLR— >0— Ti-ifl^O- 

,„ 1)( n ,. T/ ,. #V (Cl +C 3 ) 2 (C1 - C 3 ) 

(2.146) hm F = c 3 , lim — = - — - „ „. -, 

h L R^O- h LR ^o- dh LR (cf + c^)w R 

and 

(2.14c) lim *^ = lim dw ~ ^ 



Iilr^o- dh LR ci + c 3 ' h LR ^o- dh LR C1 + C3 

Remark 2.2. The limits found in (2.14a) and (2.14c) coincide with those in (2.10a) and (2.10d) 
respectively. This implies that, except at those points where V = 0, the function u- = (v-,W-) 
is of class C 1 (with Lipschitz continuous derivatives) in the whole domain {ul 6 Hi, u R 6 H3 | 

Case l)-b: Suppose that ul G Hi and u# G Hi. Here one has, by definition 

(2.15) "-^r = — \C\Wl + C\W R + f/j f^) = 1 . 

zc\ z zci 

According to our formulation (1.16) of the initiation criterion, the solution will take its values 
in phase 1 region only in case h^ R < u>m, while a phase 3 state will appear in the solution if 
h LR exceeds %. We distinguish between these two situations. 

Case l)-bl: Suppose moreover that —1 < Hlr < Wm- 

We seek the solution in the form of three constant states separated by a —c\ contact wave 
and a C\ contact wave: 



(2.16) u(t,x) 



ul for x < —c\t, 

for — c\t < x < Cit, 
u R for x > c\t. 



The intermediate state u* = (v*,w*) € Hi must satisfy the jump conditions: 

-ci(w* - w L ) + t>* - v L = 0, c 1 (w R - to*) + v R - v* = 0, 
which lead us to explicit expressions 

(2.17) t>* = v L + ci{h LR - w L ) and w* = h LR . 

Because of the assumption — 1 < h^ R < wm, it is immediate that w* belongs to H\. Note that 
no solution taken its values in H exists when h^ R < — 1. 
For further reference, we state: 

Lemma 2.3. Consider the Riemann problem (1.1), (2.1) in case l)-bl i.e. when Ul 6 Hi, 
u R € H3 and —1 < hL R < wm- Then the admissible weak solution to these problems is given 
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by (2.16), (2.17). The state u* = (v*,w*) is a C°° function of the initial states ul and ur. 
Moreover, when h^R tends to wm, one has: 

(2.18a) lim v* = v L + ci(w M - w L ), lira w * = w M , 

and 

(2.186) lim i^^(u L ,h LR ) = c l7 lim 



h LR ^w M <%LR " hr.n^wZ. dh] 



J M 



Remark 2.3. 1) It is of interest to note that the assumption —1 < Hlr < wm in Lemma 2.3 
is always fulfilled if both ul and ur belong to Hi and \ur — «l| is small enough. This is clear 
in view of (2.15). 

2) In [2], an initiation criterion was introduced in the case l)-bl. Indeed, instead of the 
solution (2.16) containing no phase boundary, the criterion in [2] selects in some cases a solution 
containing two phase boundaries (cf. (2.19) below). 

Case l)-b2: Suppose now that hLR > wm- 

According to our initiation criterion (1.16), the solution must contain (at least) one phase 3 
state. We seek the solution in the form 



(2.19) u(t,x) 



where U\,U3 £ TL\ and U2 £ Ti.3 and — cs < V < < V < c 3 . The jump conditions read: 

-c 1 (w 1 - w L ) + v x - v L = 0, ci(w R - wz) +v R -v 3 = 0, 
(2.20a) { V'(w 2 -w 1 )+v 2 -v 1 =0, V (v 2 - v x ) + c\w 2 - c\w x = 0, 

V(w R - w 3 ) + v R -v 3 = 0, V(v r - v 3 ) + c\wr - c\w z = 0. 



' ul 


for 


X < —Cit, 


Ui 


for 


- at < x < V, 


< u 2 


for 


V < x < V, 


u 3 


for 


V < x < at, 


> ur 


for 


x > at, 



They must be completed by the kinetic relation along the lines x/t = V and x/t = V: 
(2.206) 



-(fci - k 3 )(w M w m - w 1 w 2 ) = -4>(-V) 



-(fci - kz)(w M w m - w 2 wz) = 4>(V). 

It can be shown that, in fact, V' = V (cf. our previous calculations, (2.7)). As for case l)-al, 
it can be checked that (2.20) determine uniquely the admissible solution. We omit the details 
and simply state our result of continuous dependence. 

Lemma 2.4. Consider the Riemann problem (1.1), (2.1) in case l)-b2, i.e. when Uj_, € 7i\, 
ur 6 7i\ and h^R > %. Then the admissible weak solution of this problem is given by (2.19), 
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(2.20). One can consider the states m, u 2 , U3 and the speed V as functions o/ul and ur. Then 
the functions Ui and V are C 1 (with Lipschitz coninuous derivatives) and one has 

(2.21a) lim V = 0, lim vj = Vl + C\(wm — Wl) for j = 1,2 or 3, 



(2.216) 

and 

(2.21c) 



lim w\ = lim u> 3 = wm, hm w 2 = j-wm, 

hLR >W M h LR ^w M h LR ^w M K3 



lim 



ci, 



lim 



dw-: 



flLR^WM dflLR h LR ^W M dhj^Ji 



1 for j = 1 or 3. 



Note that the limits found in (2.18) and (2.21) for the functions and and Vj and wj 
(for j = 1 or 3) coincide. Hence if in case l)-bl, we set 

(2.22) Vj = u*, Wj = w* for j = 1 or 3, 

then the functions vj and Wj are globally of class C 1 with Lipschitz continuous derivatives in 
the whole range of values {ul £ Hi, £ H3}. 

From Lemmas 2.1 to 2.4, we deduce the following property of continuous dependence of the 
solution of the Riemann problem. 



Theorem 2.1. Consider the admissible weak solution to the Riemann problem (1.1), (2.1) 
described in Lemmas 2.1 to 2.4. Then the states and the wave speeds in the solution are locally 
Lipschitz continuous functions of the initial constant states ul and ur. As a consequence, if 
ui(-,0) and U2(-,0) are two Riemann initial data for system (1.1), the corresponding admissible 
solutions ui and u 2 satisfy the following L 1 continuous property: 

(2.23) / \u 2 {t,x) -ut(t,x)\dx < 0(1) / \u 2 (0, x) - «i(0, x)\ dx 

J A J A—c\t 

for all A < B and t > 0. 



We now turn to the Riemann problem in the half space ]0, oof, i.e. problem (1.1), (2.2). Two 
cases must be distinguished: 
Case 2)-a: Uo £ H3, 
Case 2)-b: u G Hi- 

Case 2)-a: Suppose that uq £ H3. 

According to condition (1.20)j, the solution must contain a phase boundary if and only if 
ho < w% where 

(2.24) h = w + — 

and is the critical value for initiation introduced in Section 1. We recall that wq > w m and 
wZ > w m . 
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Case 2) al: Suppose moreover that /io > 

Then the solution u to (1.1), (2.2) must stay entirely in phase 3, so we seek u in the form 



,.cr 



( u* for x < c 3 t, 
(2.25) u(t,x) = \ 

I no tor x > c 3 t. 

To satisfy the boundary condition (1.17a), we must have 

(2.26a) u* = 0. 

From the Rankine-Hugoniot relation c 3 (wo — w*) + v — v* = 0, we deduce w* : 
(2.266) w* 

The solution is completely determined by (2.25), (2.26). 



w -\ = h . 

c 3 



Case 2)— al: Suppose now that h < w^. 

In that case, the solution must contain a phase boundary, so we set 







f 

U- 


for 


x < Vt, 


2.27) 


u(t, x) = < 


u + 


for 


Vt < x < c 3 t, 








for 


x > c 3 t 



In view of the boundary condition (1.17a), one has 

v- = 0. 

We determine V, W- and u + = (v+,w + ) by writing the Rankine-Hugoniot relations satisfied 
along the lines x/t = V and x/t = Ci, as well as the kinetic relation along x/t = V. By 
calculations similar to those made in case l)-al, we obtain the following formulas: 



(2.28a) 



0, 



W- 



C 3 + V 

cr + V 



h , 



(2.286) v + = v + c 3 w — c 3 h = — c 3 h , w + = — h , 

c 3 -V c 3 -V c 3 -V 

where V is given by the following implicit equation: 

, OOQ ^ h-k 3[ (c 3 + V)(c 1 -V) 2 , 

(2.28c) — ^— {w M w m - {ci + v){c3 _ v) K} = <KV)- 

These formulas determine the solution in this case. 

Case 2)-b: Suppose that Uq € Hi: 

According to condition (1.20)^, we have to distinguish between two cases. Here ho = wo + — ■ 

Cl 
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Case 2) bl: Suppose moreover that ho < w^- 
Then the solution stays entirely in phase 1: 



(2.29) 

where 
(2.30) 



u(t, x) 



and 



for x < C\t, 
u for x > C\t 



Cl 



+ w = ho e Hi . 



Case 2)-b2: Suppose moreover that h > w^- 
Then the solution contains a phase boundary, i.e. 



(2.31) 



U- for a; < Vt 
u(t, x) = < u + for < x < cit 
u for x > c\t, 



where U- G Hz, u + G Hi and V e]0,ci[. The states u_, u + and the speed V are uniquely 
determined by (1.17a), the Rankine-Hugoniot relations and the kinetic relation. We omit the 
details. 

Finally, we conclude with the result of L 1 continuous dependence for the Riemann problem 
in a half space. 



Theorem 2.2. Consider the admissible weak solution of the Riemann problem (1.1), (2.2) 
described by cases 2). The states and the wave speeds in the solution are Lipschitz continuous 
functions of the initial state Uq. As a consequence, if u' and u' Q ' are two Riemann data for the 
system (1.1) in the half space, then the corresponding admissible solutions u' and u" satisfy 




u'(t,x) -u"(t,x)\dx < 0(1)(B + at - m&x(0, A - at))^ - u' ' 



for allO < A < B and t > 0. 



3. EXISTENCE VIA GLIMM'S SCHEME: STABILITY 

This section and the following one deal with the application of the random choice method, 
introduced by Glimm [15] for hyperbolic problems, to the system of mixed type (1.1). Our main 
result establishes existence of a class of admissible weak solutions to (1.1). This serves to justify 
the formulation of the Cauchy problem proposed in Section 1. 

Based on successive resolutions of Riemann problems, Glimm's method yields a sequence of 
approximate solutions for the Cauchy problem associated with (1.1). Our goal is to prove the 
convergence of these approximate solutions to an admissible weak solution to the problem, in 
case the initial data is a small BV perturbation of a single propagating phase boundary. The 
main result of this section, Theorem 3.1, yields the stability of the scheme in the BV norm. This 
guarantees its convergence in the L 1 norm to a function of bounded variation, which indeed is 
a weak solution of (1.1). Showing that this function is an admissible solution requires a more 
detailed analysis which is performed in the next section. 
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It is emphasized that a phase boundary is a wave with (necessarily) large strength. Our 
result of stability here is related to the ones obtained by Chern [5] and Schochet [37] who 
treated Glimm's scheme with large data for strictly hyperbolic systems. 

We consider the system (1.1) on the whole line (x G R) with the following Cauchy data: 



(3.1a) 



u(0, x) = uq{x) 



u° L (x) = (v° L (x), w° L (x)) for x < 0, 
^ u° R (x) = (v° R (x),w R (x)) for x>0. 



The functions u° L G BVi oc (R-,H) and u° R G BVi oc (R+,H) are assumed to be close to two given 
constant states u* L 



(3.16) 



{v* L ,w* L ) and u* R = (v* R ,w* R ) respectively, i.e. 

\ U L ~ U*l\\bV(R-) + \\U°R ~ U*r\\bV(R+) < 1- 



For definiteness, we consider the case that u* L G Hi and u* R G H3. We are assuming that 
the Riemann problem with data u* L and u* R is solved by a unique wave, i.e. a single phase 
boundary but not contact discontinuity. This assumption allows us to focus our attention on 
phase boundaries which are the main difficulty regarding system (1.1). Let u* be the solution 
of this Riemann problem; for some speed V*, one has 



(3.2a) 



u*(t, x) 



for x < V*t, 



for x > V*t. 



In the case of a characteristic phase boundary, i.e. when V* = C3, we will restrict ourselves 
to the case that no strong wave arises from a perturbation of the initial states u* L and u* R . 
According to Lemma 2.1 of Section 2 (cf. formulas (2.10a) and (2.10b)), this is the case under 
the following assumptiom: condition: 



(3.26) 



If V* = C3, then w* L = and w 



R 




Note that (3.2b) implies that h^R = and vr = vl — y/ c^c// (03) . So a Riemann problem with 
initial data in the neighborhood of ul and ur takes its values in the same neighborhood. 

We shall prove that problem (1.1), (3.1) admits an admissible weak solution which has the 
following structure: 



(3.3a) 



u(t, x) 



u L (t,x) for x<x(t), 



u R (t,x) for x>x(t), 



where 



(3.36) u L E L^([0,oo[,B^(J2,Wi)) and u R G ([0, oo[, BV(R,H 3 )) 

and 

(3.3c) X G^ 1 o ' c oo ([0,oo[, J R) and ^ G BV loc ([0, oo[,R). 
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Setting 

u*(t,x) = l U * L f ° r X<xW ' 
\ u* R for x > x(t), 

we shall also show that for all times T > 

(3.4) \\u(T) - u*(T)\\ L ^ {R>n) +TV R (u(T)-u*(T)) + TV T (^ - V*) « 1. 

The solution will be obtained as the limit of approximate solutions to (1.1), (3.1), having a 
structure similar to the one described by (3.3), (3.4). These approximate solutions are given by 
Glimm's scheme which we now describe. 



Let r > and h > be time and space mesh sizes satisfying the CFL condition tc\ < h. 
The ratio A = h/r is taken to be a constant. Let {a n } n >i be an equidistributed sequence with 
values in the interval ] — 1, 1[. We define u h (0,x) by L 2 -projection from the data u : 

(3.5a) u h (0,x) = — / uo(y) dy f° r x S [mh, (m + 2)h[ with m even. 

Note here that \u h (0, x) — u* L \ <C 1 for x < and |ti h (0, x) — <C 1 for x > 0; also u h satisfies 
u h (0, x) G Hi U H3. If u^ 1 is known up to the time t = nr — 0, we define u h {riT + 0,x) by a 
random choice projection using a n : 

(3.56) u h {riT + 0, x) = u h {nr - 0, (m + 1 + a„)/t - 0) 

for x G [m/i, (m + 2)/i[ with m + n even. Then the approximate solution u^ 1 in the strip 
{tit < t < (n + l)r} is computed by solving the Riemann problems for system (1.1) at each 
center x = mh with m + n even. 

As a consequence of our result of stability below, this construction indeed makes sense and 
yields u h (t,x) for all times t > 0. In particular, because of the assumption (3.2b), the values 
u h (t,x) stay in the neighborhoods of u* L or u* R . This implies that case l)-b2 of Section 2 will 
never occur here. The possible wave structures of the Riemann problem used in the construction 
of u h can be listed. 

Remark 3.1. As a very first step toward a general proof of convergence of u , we may consider 
the case when vP L and u R are constant, equal to u* L and u* R respectively. In that case, u h can 
be computed explicitly and consists for each time t of a single phase discontinuity connecting 
u* L at the left to u* R at the right. The position of the phase discontinuity, say ^(t), is shifted 
to the left or to the right (depending on a n and the speed V*) at each time nr. This is typical 
behavior for Glimm's scheme, which, as is well known, does not produce any numerical diffusion 
of the discontinuities. Using only the equidistributedness of {a n }, it is an easy matter to show 
that x h {t) converges to V*t = %(t) for each time t > 0. 

According to the technique of Glimm, the first step toward a proof of uniform BV stability 
for the scheme consists of studying wave interactions between Riemann wave patterns. Let 
R(u£,u r ) be the solution of the Riemann problem with initial data ut at the left and u r at the 
right (cf. Section 2 for the explicit construction). The wave strengths are defined first in case 
V* 7^ C3. We denote by £i(ug, u r ) and £2^1, u r ) the strengths of the left and right contact waves 
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in R(ue,u r ) respectively. By convention £ 2 (ui,u r ) = when R(ue,u r ) contains a supersonic 
phase boundary so that no right contact discontinuity is present. We denote by £o(ue,u r ) the 
strength of the phase boundary in R(ui, u r ) in case where ug and u r are in diffferent phases. By 
convention, the strengths are always measured in terms of the jump of the variable w across the 
wave under consideration. When a phase boundary is present, we denote its speed by V(ug, u r ). 

Consider now the case that V* =03. We define £q, £1, £2 and V in the same way as above, 
except in case that the Riemann problem R(ue,u r ) admits a supersonic phase boundary. In 
this latter case, we virtually split the phase discontinuity into two distinct waves and set: 



£ 2 (u e ,u r ) =w r - w+\ hir=Q - h er ^-\ her=0 = w r - J - 2c 3 h er /(c 1 + c 3 ), 



and 



„ . . 1 , dw + 1 / d)'{cz) ,1/ \ 

£ (u e ,u r )=w + \ her=0 + he r — \ her=0 -w- = J—— + 2c 3 h er /( Cl + cs)-w., 

where w_ and w + are the values taken by the solution of R(v,£, u r ) at the left and at the right 

of the phase boundary respectively and Ji£ r is given by (2.3). We recall that ^ + \ h( _ is given 

by Lemma 2.1. In other words, we extend the definition of £ and £2, known for h^ r > 0, to 
negative values of hg r so that their extensions are of class C . 
From the results in Lemmas 2.1 and 2.2, one easily checks that 



(3.6) 



' the functions £i,£2, £0 and V are Lipschitz continuous 

functions of their arguments (v,£,u r ), the £i 's are of class 
k C 1 (with Lipschitz continuous derivatives) away from V = 0. 



In (3.6) and all the sequel, only states which are close to either u* L or u* R are considered. 
The wave interaction estimates are derived in the following lemma: 

Lemma 3.1. Consider states ue, u p and u r which are close to either u\ or u* R . 

1) If ui 6 Hi, Up S Tii and u r 6 H3, then for j = 0, 1, 2 we have 

(3.7a) £j(ue,u r ) = £j(ue, u p ) + £j(u p , u r ) + 0(l)\£2(ue, u p )\ 
and 

(3.7b) V(u e ,u r ) = V(u p ,u r ) + 0(l)\£ 2 (ue,u p )\. 

2) If Ui G Tii, u p € H3 and u r G H3, then for j = 0, 1, 2 we have 
(3.8a) £j(u£,u r ) = £j(u£,u p ) + £j(u p ,u r ) + 0(1) D(ui,u p ,u r ) 
and 

(3.8b) V(u e ,u r ) = V(ue,u p ) + 0(1) j|£i(u p ,u r )| + |£ 2 (ii p ,u r )| j, 
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where we have set 



Si(u p ,u r )\ 



if 



V(u e ,u p ) < c 3 , 



(3.9) D(ue,u p ,u r ) 



£i(u p ,u r )\ + (V(ue,u p ) - c 3 )\£ 2 {u p ,u r )\ 



if 



V(u e ,u p ) > c 3 . 



3) if either ug, u p and u r € TL\ or ug,u p and u r £ 7is, then for j 



1,2 we have 



(3.10) 



£j(ue,u r ) =£j(ug,u p ) +£j(u p ,u r ). 



Remark 3.2. 1) Estimates in Lemma 3.1 mainly contain linear interaction terms instead of 
quadratic ones as is the case in [15]. Linear error terms were previously found useful by [5] 
and [37] to treat strictly hyperbolic systems with large data. In (3.9), the interaction term is 
proportional to the angle between the c 3 contact discontinuity and the phase boundary. Such a 
term was used by Liu [34] to analyze non-genuinely nonlinear systems of conservation laws. 

2) When V* ^ c 3 , the derivation of the estimates in Lemma 3.1 requires only the Lipschitz 
continuity of the functions £j and V, which is exactly the regularity available in general. 

3) The smallness condition (on \ug—u* L \, etc) in Lemma 3.1 is necessary only to prevent initiation 
of a new phase when solving a Riemann problem with data in a single phase. 

Proof of Lemma 3.1. We first give the proof of (3.7), and then that of (3.8). The proof of 
(3.10) is trivial. 

In view of the results of Section 2, the functions £q, £\, £2 and V are (at least) Lipschitz 
continuous functions of their arguments. Hence the formulas (3.7) will follow easily if we check 
that 

(3.11a) £j(ug,u r ) = £j(ug,Up) + £j(u p ,u r ) 



hold whenever £i{ug, u p ) = 0, i.e. when there is no right wave in the left wave packet R(ug, u p ). 
But this last statement is obvious because the left-waves in R(ug,u p ) and R(u p ,u r ) are associ- 
ated with a linearly degenerate characteristic field. Such waves can be superimposed without 
any interaction and the wave strengths are simply summed up, cf. (3.11a). The speed of the 
phase boundary remains unchanged, cf. (3.11b). (These facts can be checked directly from the 
analytical expressions in Section 2.) The proof of (3.7) is completed. 
We now prove (3.8). We notice first that (3.11a) as well as 

(3.11c) V(u e ,u r ) = V(ug,u p ) 

do hold provided that D given by (3.9) vanishes. 

Namely, when V(ug,u p ) < c 3 and if D(ui,u p ,u r ) = 0, that means that the right wave 
packet does not contain a — c 3 contact wave. In that situation, the two wave patterns can be 
superimposed, without any interaction. When V(ug,u p ) > c 3 and if D(ug,u p ,u r ) = 0, then 
the right wave packet does not contain a — c 3 contact wave and, moreover, either it also has no 



and 



(3.116) 



V(u e ,u r ) = V(u p ,u r ) 
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c 3 contact wave or the speed V(ue,u p ) equals C3. In both situations, the left and right wave 
patterns can be superimposed. Again, there is not interaction. This proves (3.11a) and (3.11c). 

When V(u£, u p ) < C3, estimates (3.8) follow from (3.11a), (3.11c) and the Lipschitz continuity 
of £j and V. When V(ue,u p ) > C3, since V is away from V = and according to the results 
in Section 2, the functions £j and V are of class W 2, °° . This allows us to apply the classical 
lemma of division (e.g. [Ho]) and again to deduce (3.8) from (3.11a), (3.11c). The proof of (3.8) 
is completed. □ 

We will use the technique of Glimm to deduce from Lemma 3.1 the result of BV stability of 
the scheme. We refer to [15] and [16] for the terminology we use here. At this stage, we have 
to define functionals to control the total variation of the solutions. The choice we propose here 
is motivated by the form of the terms of interaction found in Lemma 3.1. Note that the phase 
boundary which is a "strong wave" are treated separately from the "small waves" . 

The (t, x)-plane is divided into a set of diamonds A mj „ with centers (nr, mh) (n + m even) 
and with vertices 



Given a diamond A mn , we denote by un , «b, uw and us the values taken by u at the vertices 
N, E, W and S respectively. 

We give now the definition of the approximate phase boundary in u h that we denote by 
X h - R+ — ► R- First of all, it is a simple (but useful) observation that the phase boundary in u h 
is actually located at a single space position for each time t = nr. In other words, there is no 
spreading of the phase boundary. Let \ h De the piecewise linear curve which is discontinuous 
at each t = nr and coincides with the phase boundary in u h inside each slab [nr, (n + l)r[. Let 
V be the set of all diamonds that are crossed out by the phase boundary \ h ■ 

We then introduce several functionals defined on space-like curves, say J, passing through 
vertices of diamonds. Define 



where £q is the strength of the phase boundary when crossing the curve J. The functional L(J) 
bounds the total variation of u h along the curve J on both sides of the phase boundary. B( J) 
measures the jump of u h across the phase boundary. Next we define the potential interaction 
Q(A) in a diamond A by 



N = ((n + 1)t, (m + a n+ i)h), E = (nr, (m + 1 + a n )h), 
W = (nr, (m — 1 + a n )h), S = ( (n — l)r, (m + a n _i 



)h). 



(3.12a) 




the summation being on all small waves crossing the curve J, and 



(3.126) 



B(J) = \£ 



\£ 2 {u w ,u s )\ 
(3.13) Q(A) = { \£i(u s ,u E )\ 



( 



if 



if 



A = A m; „ G V and mh < x h { n ^) 



+6(V{u w ,u s ) - c 3 )\£ 2 (us,u E )\ 



V 



if 



A 



A m ,n £ ^ and mh > x ?l (nr), 



where 6: R — ► R is the function defined by 



9{y) = for y < 0, y for y > 0. 
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Finally the potential wave interaction Q(J) associated with a curve J is 
(3.14) Q(J)= 1^1+ E {\£i\ + HV-c 3 )\£ 2 \\, 

waves at the waves at the ^ ' 

left of x h right of x h 

where V is the speed of the phase boundary when crossing J and the summation being on all 
waves crossing J. Note that Q(J) is a linear functional in terms of wave strengths. 

Lemma 3.2. Let K be a sufficiently large constant. Let Ji and J 2 be two space-like curves, J 2 
being a successor of Ji . Then one has 

(3.15a) L(J 2 ) + KQ(J 2 ) < L(Ji) + XQ(Ji) 

and 

(3.156) B(J 2 ) + KQ{J 2 ) < B{J 2 ) + KQ{J 2 ). 



Proof. We need only prove (3.15) when J 2 is an immediate successor of Ji; the general case 
follows by induction. We check first the formula 

(3.16) q(j 2 )_q(j 1 )<_Iq(a), 

where A is the diamond limited by Ji and J 2 . If A 6 V and the right wave packet contains the 
phase discontinuity, then in view of (3.14), (3.13): 

Q(J 2 ) - Q(Ji) = -\S 2 (u w ,u w )\ = -Q(A) < ~Q(A). 
If A G V and the left wave packet contains the phase discontinuity, then in view of (3.14) 

Q(J 2 ) - Q(Ji) = -\£i(u s ,u E )\ - 0(V(u w ,u s ) - c 3 )\£ 2 (u s ,u E )\ 

+ le(V(u w ,u E )-c 3 )-e{V{u w ,u s )-c 3 )\ E \ e *\ 

^ * waves on the right 

side of A 

Since 9 is Lipschitz continuous and using definition (3.13) and (3.8b), it follows that 
Q(J 2 )-Q(Ji) = -Q(A) + 0(l)Q(A) £ \£ 2 \ 

waves on the right 
side of A 

= Q(A)j- 1 + 0(1) L(J!)}<~Q(A) 

where in the last inequality we have assumed that 0(1) L(J X ) < \. The condition L(Ji) ^ 1 
indeed is ensured by induction if the initial total variation is small enough. Let us content 
ourselves with checking that L(Jq) <^ 1 where Jo is the curve connecting centers of diamonds 
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on the lines t = and t = t\ = T. Namely, using the definition of the wave strengths and (3.2b), 
one gets 

L(j ) = o(i)|ry_ oo «) + tv +oo (u° r ) + 1|< - <|U~ (JM + K - u^||^ (R+) |. 

The right hand side of the above formula is small in view of (3.1b). 
If now A g" V, then from (3.10) one has trivially 

Q( J 2 ) - Q(Ji) < and Q(A) = 0. 

Henceforth, the proof of (3.16) is completed. 

We next consider L(J 2 ) - L(Ji). If A G V, then by (3.7a), (3.8a), (3.12a), (3.13): 

L(J 2 ) - L(J X ) = \£i(u w ,u E )\ + \£2(u w , u E )\ 

- \£i(u w ,u s )\ - \£ 2 {u w ,u s )\ - \£i(u s ,u E )\ - \£2{u s ,u E )\ 
= 0(1)Q(A). 

If A P, one has 

L(J 2 ) = L(Ji) and Q(A) = 0. 

This proves the formula 

(3.17) L(J 2 ) = L(J 1 ) + 0(1)Q(A). 

From (3.16) and (3.17), we easily deduce (3.15a) provided that the constant K in (3.15a) is 
large enough. 

Finally, it can be proved similarly that 

(3.18) B(J 2 ) = B(J 1 ) + 0(1)Q(A) 

which implies (3.15b) in view of (3.16). □ 

Lemma 3.2 provides a uniform bound for the total variation of u h (t) at times t = t n . Since 

TV{u h {t)) < 0(1)L(J), for all times t € [t n ,t n+1 [, 

where J is the curve lying between the lines t = t n and t = t n+1 , we obtain a uniform control 
of the total variation of u h for all times. Let us define the function u h : R + x R — > H by: 



(3.19) u h {t,x) 



u* L if x < X h (t), 
u* R if x > X h (t), 



where \ h 18 the approximate phase boundary associated with u h . From Lemma 3.2, one deduces 
the following result of stability. 

Theorem 3.1. The functions u h given by Glimm's scheme applied to the mixed system (1.1) 
and the data (3.1), (3.2) satisfy the following stability estimates: 

(3.20a) TV±™(u h - u h ){t) < 0(1)N U 
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(3.206) 

and 

(3.20c) 



\u h (t) -u h (t)h-(R,H) < 0(1)(N 1 +N 2 ) 



\u h (t) - u h (t')\\ LHR , n) < 0(1)N! (\t -t'\+h 



for all times t > and t' > 0, with Ni and N 2 defined by 

{ TVL^K) + TV °°(u° R ) if V* ± c 3 , 



(3.21o) 

and 
(3.216) 



k Ty_° oo (n0 ) + TV£°(u%) + N t if V* = c 3 



N 2 = \\u° L - u* L \\ L ~ {R _ iHl) + \\u° R - u* R \\ L ^ {R+tH3) . 



By Helly's theorem, the estimates (3.20) imply that (a subsequence of) {u h } converges in 
L\ oc strongly to a function u as h — > 0. This function has bounded variation in space and 
satisfies the same bounds as u h in (3.20). It is a classical matter (Glimm [15], Liu [33], [34]) to 
check that u indeed is a weak solution to the system of conservation laws (1.1). It also satisfies 
the entropy inequality as well as the initial condition. It remains to show that u is admissible, 
i.e. satisfies the kinetic relation (cf. Section 4). 

Remark 3.3. If condition (3.2b) is violated, then perturbating of a characteristic phase bound- 
ary produces a c 3 -contact wave with strong strength. Then, one would have to deal with inter- 
action between two strong waves traveling with arbitrary close speeds. Initiation of new phases 
is possible. It is not clear whether the total variation of u h would remain uniformly bounded in 
that case. 
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4. EXISTENCE VIA GLIMM'S SCHEME: ADMISSIBILITY 

In this section, we prove that the weak solution u = limu h found in Section 3 by Glimm's 
scheme does satisfy the kinetic relation (1.11). This establishes that u is an admissible weak 
solution to our problem and leads to the desired result of existence and stability. 

First of all, we notice that the kinetic relation (1.11) is formulated in a pointwise sense, 
more precisely (1.11) must hold almost everywhere with respect to the Hausdorff measure Hi. 
However from the results in Section 3, we only have that u h converges to u at almost every point 
with respect to the Lebesgue measure on R + x R. This latter property is thus not sufficient to 
pass to the limit in the kinetic relation. 

We prove in this section a result of pointwise convergence for the approximate phase boundary 
X , cf. Theorem 4.1. This result is derived by using the techniques introduced by Glimm and 
Lax in [16]. The focus of [16] was the case of a strictly hyperbolic system of two conservation 
laws with small data. Extensions of the results in [16] can also be found in the papers of 
Di Perna [10] and Liu [34]. In our situation, we have a (special case of) a system of mixed type 
with large data. 

Next in Theorem 4.2, we prove that the above result is sufficient for the passage to the limit 
in the kinetic relation, assuming that the speed of the phase boundary V* does not vanish. 

Let us consider the phase boundary x h '- R+ — > -R in the approximate solution u h . The 
function x h is discontinuous and piecewise linear. It jumps up to a distance of ±2h at each 
time step. It is easy to verify the following lemma. 

Lemma 4.1. The function x h - R+ — ► R satisfies the following uniform estimate: 
(4.1) |^(t)_^( t ')|< h t -t\+2h for 0<t<f. 



By Ascoli's theorem, the sequence {x h } must converge on each compact set in the uniform 
topology to a function \ G W l( £°{]$),oo[,R). The next lemma gives a bound for the total 
variation of the functions \ h '■ R+ ~ > R defined by 

dx h (t) 

(4.2) x (t) = — ^ — (constant) on each interval [nr, (n + l)r[. 

From an analysis of the waves crossing the phase boundary, we prove the following result. (The 
proof is given after the statement of Theorem 4.1.) 

Lemma 4.2. For all times T > 0, one has the uniform estimate 

(4.3) TV T ( X h ) < O(l){TV^ T/x) _ 2h (vl-c lW l)+TV^ /X)+2h (u R ) + N}, 

where TV = if V* ^ c 3 and TV = \\u° R - «2lU°°(-(r/A)-2h) + \W R - u* R \\ L ^ {{hT / x+2h) if V* = c 3 . 

Hence, from Lemmas 4.1 and 4.2, the equidistributedness of the sequence {a n } and arguments 
in [16], we deduce the following pointwise convergence property. 

dx 

Theorem 4.1. The functions x h and x h converge to the functions x and —— respectively in the 

dt 

following sense: 

(4.4) \\x — X h \\L°°(]o,T[,R) — ► when h — > 0, for all T>0 



31 



and 

(4.5) x h (t) -» for all times t€R+\E, 

where E C i?+ is an at most countable set. 

We give first the proof of Lemma 4.2 and then the one of Theorem 4.1. 

Proof of Lemma 4.2. Let T be fixed and let N be such that Nt <T < (N + 1)t. Recall that 
P is the set of all diamonds which contain a part of the phase boundary. Let J be the spacedike 
curve which limits the domain of dependence of the diamonds in T> with centers below the line 
t = Nr. For each time t = nr, n = 0, 1, . . . ,N, J encloses a finite number of diamonds that we 
denote by A™ for m = 1,2,... , N + 1 — n. They are ordered increasingly. We define m{n) to 
be such that A™ (n) £ V. 

By Lemma 3.1, the speed of \ h at the time nr is estimated from the its value at time (n — l)r: 

(4.6) x(nr + 0) = ± h { (n - l)r + 0) + 0(l)|£(A n , m(n) )|, 

where |£(A n m ( n ))| represent the strength of the waves entering the diamond. By summation 
with respect to n = 1, . . . , N, we obtain 

JV 

r(AT+l)r, 



TV^+^Oc*) = £ IxVr + 0) - *( (n - l)r + 0)| 

n=l 

JV-1 

= 0(l)^|£(A nim(n) )|, 



n=0 

which bounded by the total variation of u h measured along both sides of the phase boundary. 
By using conservation laws for wave stregnths, as made in [16], one could check that the total 
variation of u h along this curve is bounded by the initial total variation. Thus we have 

N+l 

(4-7) TV {N+1)T (x h ) = O(l)Y J \£(^0, m )\, 

m=l 

with 

m(0) 

(4.8a) £ |f (A , m )| = 0(1) TV° (T/x) _ 2h (v° L - c lW ° L ) 



m=l 



and 



,, s ., ^ .... .. ro(l)TVf™«V H ) if y. #C3 



Combining (4.7) and (4.8) gives (4.3). The proof of the lemma is complete. □ 
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Proof of Theorem 4.1. In view of Lemma 4.1 and Ascoli's theorem, we have the convergence 
result (4.4). In view of Lemma 4.2, the total variation of x h on a compact set [0, T] is uniformly 
bounded. By extracting again a subsequence, Helly's Theorem gives 

(4.9a) x h {t)^k(t) for all times t > 

and 

(4.96) TV^{± h ) -» i(t) for all times t > 0, 

where k e SVJ oc ([0, oo[,J2) and £ € L^.([0,oo[,J2 + ). We define E C i?+ as the set of all 
points of discontinuity of the function £. This set is at most countable since the function £ is 
non-decreasing (and so has bounded variation). 

We are going to prove that (4.4) holds with the above choice of set E. Let t be in R + \ E, 
and let e > be so small that 

(4.10) TV?+/(x fc ) < e- 

This is possible because t ^ E. Then, in view of (4.9a) and (4.10), we have 

(4.11) k(t) - e < X h {s) < k(t) + e for s e]t - e, t + e[. 

On the other hand, we know that the curve x h has the slope x h {t) on the interval [nr, (n+l)r[3 t 
and jumps by ±2h at times (n + l)r. The slope x h of x h is "controlled" by inequalities (4.11), 
while the jumps of are determined by the given sequence {a n }. 

Let n 1 and n" be two integers such that (n' — l)r < t' < n'r and n"r < t" < (n" + l)r, 
where t — e < t' < t" < t + e are given. We set 

0+ = {m/m integer, ri < m < n" and a m < (k(t) — 

and 

7- 

rT = integers, n < m < n" and a m > (fe(t) + 

In view of (4.11) and between times n'r and n"r, the curve x' 1 has at least #ri* jumps to the 
right and at most n" — n' — to the left, thus we have 

(4.12a) X h i.n"r) - x /l (n'r) > (2#ft* - n' + n")h. 

Similarly for Q* we get 

(4.126) X h (n"r) - x h ip!r) < (2# - ri + ri')h. 

But the equidistributedness of {a n } means that 

when h — > 0. 
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Combining (4.12) and (4.13) and letting h — > yield the inequalities 



x(t") 



x(0 > (*:(*) - e) (* 



// 



and 





Letting e go to zero in (4.14) yields 



dx{t) 



k(t). 



dt 



The proof is complete. □ 

Remark 4.1. Estimate (4.3) of Lemma 4.2 makes clear that only the c\ waves located at the 
left of the initial phase discontinuity and the ±03 waves located at the right of the initial phase 
discontinuity contribute to the change in speed of the phase boundary. 

We finally prove that the result in Theorem 4.1 is sufficient for the passage to the limit in 
the kinetic relation. 

Theorem 4.2. Suppose that V* ^ 0. Then the limit function u given by Glimm's scheme 
satisfies 



Hi-almost everywhere on the set B su b(u). 

Equality (4.15) is understood as equality between Borel measures on R + x R. Here B su b{u) 
(according to the definition of Section 1) is the set of all points of approximate jump of u 
associated with subsonic phase discontinuities. In view of the formula of Section 1, it is clear 
that, when V* ^ 0, (4.15) is equivalent to the formulation (1.11) of the kinetic relation. The 
case V* = could in principle be treated by the same technique but this would require further 
analysis. 

Remark 4.2. 1) The pointwise convergence property of Glimm's scheme was already used in 
LeFloch-Liu [30] to derive an existence result for nonlinear hyperbolic systems in nonconserva- 
tive form. 

2) If V* = 0, x ma y vanish and then relation (4.15) is not sufficient to uniquely characterize 
the solution (e.g. of the Riemann problem). 

Proof of Theorem 4.2. For all times t > 0, we introduce an approximate normal u h (t) 



(4.15) 



d t U{u) + d x F(u) = -v t 4>{ )£c= x (t) 



X 



(!/*(*),!/£(*) )ty 



v^{tf + v h x {tf = l, X \t) 



and > 0. 
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Similarly, from x(t), we define u(t) = (z^(i), v x (t)). According to the notation of Section 1, we 
have in fact u(t) = v(t,x(t) )• First of all, we claim that 

(4.16) u^( x h )S x=xh - M(^)5 X = X 

in the weak-star topology of bounded Borel measures on R + x R. 

Since x h satisfies (4.3) and the right hand side of (4.3) is small by the assumption (3.1b), 
the function x h has small total variation. When V* ^ 0, we can ensure that x h is bounded 
away from zero uniformly with respect to h. In view of properties (1.8), the function 4> is (at 
least) continuous in the range of values taken by x h ■ This fact combined with the result of 
convergence (4.5) gives 

(4.17) ^(t)0(x h W)->^(*)^(^(*)) ^ all t€R+\E, 

where E is an at most countable set. From (4.17) and the uniform convergence of x h to x 
Cf. (4.4), we deduce that 



if (t) 4>{x h (t) ) 9{ X h (t) dt^j v t {t) 4>(^(t) ) 9{ X (t) 



dt 



for each continuous function 6: R R with compact support. This proves (4.16). 
By construction, the approximate solutions u h satisfy the kinetic relation 



(4.18) 8 t U(u h ) + d x F{u h ) = -v$<l>(— L) 5 x=x u 

ifi-almost everywhere on the set B su b(u h ). We claim that using (4.16) we can pass to the limit 
in (4.18) and obtain 

(4.19) d t U(u) + 8 x F{u) = -u t <f>(-^-) 5 X=X 

v x 

^-almost everywhere on the set B su b(u). 

The left hand sides of (4.18) and (4.19) are treated easily since they have a (divergence-like) 
conservation form. In particular, we have 

(4.20) d t U(u h ) + 3 x F(u h ) -> d t U(u) + d x F{u) 

in the weak star topology of bounded Borel measures on R + x R. 

In case V* < C3, (4.18) is satisfied on the whole space R + x R and so the desired result (4.19) 
is an immediate consequence of (4.18), (4.16) and (4.20). 

When V* > C3, nothing has to be proved since no kinetic relation is imposed then. 

The final case V* = C3 is treated as follows. We note that one can find two Lipschitz 
continuous functions (j){V) and (f)+(V) defined for V in a neighborhood of C3 such that the 
kinetic relation (e.g. for u h ) is equivalent to the two inequalities 

(4.21a) d t U{u h ) + d x F(u h ) < -^ + {-\) 5 x=%h 



35 



and 

(4.216) d t U(u h ) + d x F(u h ) > -i%4>_{-\) 5 x=xh , 

where (f>± are chosen in such a way that 

4>+(V) = 4>-{V) = <j>(V) for V<c 3 

and 

4> + and 4>- are Lipschitz continuous with: 4>-{V) < 0+(V). 

Namely, this is possible since (4.21a), (4.21b) when V < c 3 give back the kinetic relation; 
while for V > c 3 (4.21a), (4.21b) are trivially satisfied provided that the entropy dissipation 
in the supersonic case remains in the interval [4>-(V),(fi+(V)]. In this latter case, the entropy 
dissipation across the phase boundary, say 4>(V), is the following (cf. the notation of Section 2): 

4>{V) = \{k\- k 3 )(wMW m - WRW-) = \{ki- k 3 )(w M W m - ° 3 J_ Yf WRhLR), 

where V = V^Hlr) > C3 is a root of the equation (2.12b). By (1.8b) and (1.9b), we have 
Jim 4>(V) = -(fci - k 3 )w M w m = tp(c 3 ) = Jim 4>{V). 

This proves the continuity of the entropy dissipation at V = c 3 . Moreover <p is clearly Lipschitz 
continuous in view of Lemma 2.2. 

Hence, for <p±(V) suitably chosen and V — c 3 sufficiently small, the entropy dissipation 4>(V) 
remains in the interval [(j)-(V),<p + (V)]. 

It is clear that (4.16) still holds if cf> is replaced by 0_ or </> + , i.e. we have in the weak-star 
topology: 

(4.22) ^(^)W^*^(§)**=x- 
Then (4.20) and (4.22) used in (4.21) yield: 

d t U(u) +d x F(u) < -u t 4> + (-—)5 x=x h 

v x 

and 

d t U{u) + d x F(u) > -v t 4>_( — -)5 x=x h 

which give (4.15). The proof is complete. □ 

We summarize in the following theorem the results obtained along Section 3 and in the 
present section. 

Theorem 4.3. Consider the mixed system (1.1) with an initial condition which is a small 
perturbation in the BV norm of a single propagating phase boundary with speed V* . Suppose 
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that V* 7^ and condition (3.2b) is satisfied ifV* = C3. Then Glimm's scheme for this problem 
converges to an admissible weak solution which has the structure described in (3.3), (3.4). 

Remark 4.3. 1) Note that the proof of Theorem 4.2 uses the property that the entropy 
dissipation across a contact discontinuity is identically zero. 

2) We believe that Theorem 4.3 could be extended to a finite number of propagating phase 
boundaries. Also the restriction V* ^ is only a technical assumption and could be removed 
by using other techniques from [16]. 

3) However, there is a main obstacle to a general result of existence of BV solutions for (1.1). 
Indeed, for arbitrary large data, the phenomenon of initiation of new phases arises, and it is an 
open problem to derive a uniform bound on the total variation in that case. 
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